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For a number of nanophotonic applications, the conditions under which nanoparticle systems 
have polarization-invariant optical responses remains a largely unresolved issue. We study the 
optical responses of structures that have an n-fold (n > 3) rotational symmetry with respect to 
an axis parallel to the propagation direction of the incident plane wave. Such symmetric systems 
are considered as constructed from small {R <^ A) spherical particles that can be approximated 
by electric and magnetic dipoles, which allows us to employ the coupled dipole approximation 
method. It is rigorously proven that all the cross-sections (extinction, scattering and absorption) 
are independent of the angle of polarization of the linearly-polarized plane wave. It is shown that 
such polarization-independent responses are purely due to the rotational symmetry of the whole 
structure and that it is a generic geometric characteristic of any system possessing such symmetry. 
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I. INTRODUCTION 

The current surging interest in various applications 
of nanoscale light-matter interactions, including biosens- 
ing [1-3], nanoantennas [4, 5], photovoltaic devices [6] 
and many others, has triggered enormous effort into the 
old and fundamental problem of manipulation of a par- 
ticle's scattering and absorption characteristics [7, 8]. 
In the recently emerging fields of nanophotonics, vari- 
ous novel phenomena have been demonstrated involving 
interaction of nanoparticles with light, such as super- 
scattering [9, 10], control of the direction of the scattered 
light by metasurfaces [11, 12], coherent perfect absorp- 
tion of light by surface plasmons [13], Fano resonances in 
nanoscale structures [14, 15] and plasmonic oligomers [16- 
22]. At the same time, the interest in artificial magnetic 
responses that was fostered by the field of metamaterials 
has lead to the observation of artificial magnetic modes 
in nanoparticles and, since then, many related novel scat- 
tering features based on the interplay of both electric and 
magnetic responses have been demonstrated [23-29]. 

To make further breakthroughs in different applica- 
tions based on the particles scattering, there is a funda- 
mental challenge to overcome: polarization dependence. 
The dependence of an optical response on polarization 
comes from the fact that most structures have domi- 
nant ly electric responses, which are highly dependent on 
the polarization of the incident field. The simplest struc- 
ture that does not exhibit polarization-dependent scat- 
tering properties is a single spherical particle. According 
to the Mie theory the total extinction, scattering and 
absorption cross-sections do not depend on the incident 
polarization angle, although the scattering diagram will 
exhibit some angle-dependent properties [7]. However, it 
is possible to achieve polarization-independent scattering 
diagram by overlapping of an electric and magnetic dipole 
responses of a single spherical nanoparticle [28, 29], but 
such effects can be only achieved by rigorous structure 
engineering and can happen only in specific spectral 
regime. It has also been experimentally observed that 



some plasmonic oligomer structures exhibit polarization- 
independent extinction cross-sections [19, 30]. This all 
leads to the question of what the necessary conditions are 
for an arbitrary system to have polarization-independent 
scattering properties. 

Inspired by the concepts of symmetry induced degen- 
erate states in quantum mechanics [31] and modes de- 
generacy in uniform waveguides [32-34] , there have been 
some studies about the symmetry-induced polarization- 
independent scattering by clusters of particles [16, 35, 
36]. However, as far as we know, there are no rigorous 
and systematic investigations of this topic. Addition- 
ally, in previous studies, usually only the dependence of 
extinction or scattering spectra on polarization is inves- 
tigated and the intrinsic loss spectra is neglected, which 
could be quite important in its own right, particularly 
for photovoltaic devices and biological applications. 

In this paper we perform a systematic investigation on 
the optical responses of structures with an n-fold (n > 3) 
rotational symmetry, where the n-fold symmetry axis is 
parallel to the direction of propagation of the incident 
plane wave [see Fig. 1]. Such n-fold symmetry implies 
that the optical properties of the system will be identi- 
cal when rotating the whole structure by 27r/n radians. 
But, as we analytically prove, for structures with n > 3 
symmetries the extinction, scattering and even absorp- 
tion cross-sections are all identical for rotation on any 
angle. Such structures can therefore be considered as 
being polarization-independent. 

The extinction and scattering cross-sections are de- 
fined in the far-field, but absorption can be calculated 
by two independent ways - as an energy balance between 
the far field scattered and incident fields, and as inte- 
gration of losses in the near field. Both approaches pro- 
duce the same result. In the near field, the full pro- 
file of the electromagnetic field should be taken into ac- 
count, while in the far field only the leading order will 
survive. The polarization-independent absorption is then 
quite counter-intuitive because the near field profile of the 
electromagnetic field does depend on the incident polar- 
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ization yet the overall absorption cross-section does not. 
It implies that the variation of the near field with the 
incident polarization does not affect the overall integral 
absorption cross-section. Thus, the near field distribu- 
tion still follows some symmetry properties of the entire 
structure. We are able to trace such peculiar relations us- 
ing the dipole approximation because it is able to show 
the key aspects of the observed phenomena. 

The study in this paper is on symmetric systems con- 
sisting of small {R <C A) spherical particles, which can 
be approximated as electric and magnetic dipoles using 
coupled dipole approximation [37] . Given that any struc- 
ture can ultimately be decomposed into small spherical 
particles and subsequently a composite of electric and 
magnetic dipoles, the derived polarization-independent 
response should be a typical feature of any system ex- 
hibiting such an n-fold rotational symmetry. As this is 
a purely geometric feature and shows no dependence on 
the wavelength, the optical properties of the constituent 
materials or the resonances excited, we expect that it 
could be applied in various applications including sens- 
ing, imaging, photovoltaic devices and other biological 
and medical research. 



II. THEORETICAL MODEL 

A. Discrete dipole approximation 

The scattering and absorption of light by a small spher- 
ical particle {R <C A) can be approximated by that of an 
electric and magnetic dipoles with the effective polariz- 
abilities [37] 



m = 



(1) 



where p, m are the effective electric and magnetic dipole 
moments, E, H are the electric and magnetic fields act- 
ing on the particle and a^;, an are the scalar effective 
polarizabilities 



The free space dyadic Green's functions are given by 

1 

V X G%ri, r,) = = dij (uji ® n.^ - J^^)) ^ (4) 

in these we have expressed the operator uji x /^^^ explic- 
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where rji = |r — r_^-| and uji = ^ (r^ - r_^-). 

In any system constructed from small, spherical par- 
ticles the fields acting on the ^th particle will be the sum 
of both the externally- applied incident fields and the scat- 
tered fields from all the other particles. Moreover, by us- 
ing the expression for scattering from any single particle 
as given in (3), the expressions for the dipole moments 
of the zth particle are 



Mo .(3) 
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I N 1 3 
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where E^^ = E^(ri), H^^ = H°(ri) are the externally- 
applied electric and magnetic fields at r^. 



Matrix Notation 



aE = Ginai/k^ 
an = QiTrbi/k^ 



(2) 



as defined here in terms of the Mie dipole scattering co- 
efficients ai, hi. 

The scattered fields from a single particle in free 
space can be described in this approximation with dyadic 
Green's functions 



E"(r) = k^ 
H^(r) = e 



eo V ^0 ^ ^ 

G°(r,ro)m+ (v x G°(r,ro)) p 



(3) 



To describe an entire dipole system we define a state 
consisting of all dipole moments, and a state consist- 
ing the incident fields acting on each dipole, |e) 



\p) = {pi 
\e) = ( E«i 
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where the magnitudes of the magnetic fields and dipole 
moments in these states have been scaled to be that of 
the associated electric fields and dipole moments using 

scaling factors of and y/eofiQ. 

The dipole equations in (6) show that each incident 
field vector can be equated to some linear combination 
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of the dipole moments and, as such, there will exist some 
6N X 6N matrix, M(^^\ to relate |e) to \p) as 



(8) 



where a^^^^ is a diagonal matrix containing the electric 
and magnetic scalar polarizabilities, aE and an^ of each 
particle. 

The matrix in (8) is referred to as the interaction ma- 
trix. To generate a general expression for this matrix, we 
define four smaller matrices; two that couple electric or 
magnetic dipoles together, M^e^^ and M^j^\ and two 
that couple electric to magnetic dipoles and vice versa, 
and M^^^ . These matrices are constructed from 
the dyadic Green's functions as to match the concatena- 
tion of vectors in \p) and |e) [see eq.(7)] 



unchanged. This sort of symmetry is also referred to as 
cyclic symmetry. 

In the coupled dipole equation (8) any rotational sym- 
metry must be contained within the set of dyadic Green's 
functions (G)^- and g^j ) that make up the interaction 
matrix as no other terms contain information on the 
structure of the system. As such we are able to con- 
sider only the dyadic Green's functions when applying a 
general rotation operator, R^^\ to the position vector of 
every particle. Moreover, it is only the unit-dyad terms 



m G\f and gl]' that are affected because the coefficients, 
aij , bij and dij , are all independent of the position vector 
directions [see eqs.(4) and (5)]. Explicitly, it follows from 

^ (3) (3) 

the definitions in (4) that both and g\j will trans- 
form as a change of basis under any rotation applied to 
the system's position vectors 
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where aE,H is a diagonal matrix containing the electric 
or magnetic scalar polarizability of each particle and 
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The interaction matrix that satisfies (8) is then con- 
structed as 



G 



(3) 



^(3) 
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(11) 



(9) In order for (11) to follow for the case of g^j\ we must 

(3) 

also acknowledge that g^j^ can always be expanded as a 
Taylor series. This is true because the magnitude of every 

component in ^n^^ (8) uji — 1^ will be less than or equal to 

1 

one and therefore the Taylor series for (jiji Uji — ^ , 
and hence gl^\ will always converge. 

So, given that (11) holds for all G^^^ and glj \ the entire 
interaction matrix will transform in an equivalent manner 



where R^^^") is defined such that it applies R^^'' to the 
position vector of every particle in M^^^^ 



^(6iV) ^ ^(3)j(6iV)_ 
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However, for the case of a symmetric rotation about the 
^^^^ n-fold axis, Rn \ it also follows that 



III. POLARIZATION INVARIANCE AND 
SYMMETRY 

A. Rotations and Cyclic Symmetry 



rk,i 



(14) 



for each i, j G {1, ...TV} and A:, / G {1, ...TV}. Subsequently 
the dyadic Green's functions will also transform as 



In this section we are going to consider particle systems 
that have a rotational symmetry described by an n-fold 
axis. An n-fold axis is an axis in our system about 
which any number of ^ rotations will leave the system 



p(3)A(3) 



ki 

(3) 
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Therefore there raust exist some permutation matrix, 
fl.n^\ for each ^^f^such that 

We can consider ah non-zero components of H^f^^as /^^^ 
matrices and therefore nl^^^wih commute with Rn^^- 
From (16), it then fohows that Rn^'^fl-n^'^ wih commute 
with the interaction matrix 



(17) 



If we also only consider incident fields that do not break 
the symmetry of the system, that is; fields such that both 
= Bk and = hj^ for all i. A: G {1, ...N} where Rn^Yi = 
Fjt, then 



ni'^^) le) = le) 



(18) 



The dipole moments of a symmetry-rotated incident field 
can then be defined as a result of both (17) and (18). 
Specifically, for some |e), \p) such that eoQ;(*'^) |e) = 

eoa^'^^ [r^^^ |e)) = M^^^^ (^i'^^fl^f^' \p)) ■ (19) 

A particular subset of incident fields that will obey 
the result in (19) are the linearly-polarized plane waves 
travelling parallel to the n-fold axis of our system. To 
construct a general linearly-polarized wave we consider 
some orthonormal basis of two plane waves that are in- 
phase, \eh) and |e^). Then any linearly-polarized incident 
field can be expressed as 



|e) = a|e/,) +/3|e^) ; (a, /3 G R) . 



(20) 



However, if the n-fold axis of our system satisfies n > 2, 
we can choose Rn^ to be a rotation of ^ about the n-fold 
axis and then express the incident field in an alternate 

basis of ||e/^) ^Rn^^ where 

|e.) = a \eh) + hR^^''^ |e^) ; (a, 6 G R) . (21) 

The prescribed orthonormality of the original basis 
states, \eh) and |e^), requires that 
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{ey\ey) = o? -\- 2abcos{ — ) -h ti^ 
(e^le/i) = a -h 6cos( — ) = 0. 



(22) 
(23) 



In this basis the general form of any, linearly-polarized, 
incident field becomes 



|e) = (a + Pa) \eh) + ^56^1^^) \eh) 



(24) 



and we can use the relation in (19) to express the asso- 
ciated dipole moments as 



\p) = {a + M \Ph) + /36^(f^)n(f^) \ph) ■ (25) 



B. Total Scattering and Absorption Invariance 

The extinction, absorption and scattering cross sec- 
tions for a system with dipole moments \p) and incident 
fields |e) are 



eo|EO| 

k 



rlm{(e|p)} 



(26) 



eo2|E' 



(Js = (Je — CTn 



The rotational symmetry we introduced in the previous 
section is contained entirely within the inner product 
terms of the extinction and absorption cross sections, 
being (e|p), {p\a~^\p) and {p\p). We will consider these 
inner products separately. 

The first of the inner products, (e|p), can be expanded 
in terms of {eh\ and \ph) using (24) and (25) 



{e\p) = {eh\ + 2a/3a + (3^a^ + pH^fll^''^ 



The permutation matrices are then absorbed by {eh\ us- 
ing (18) 



(e|p) = {e,\a' + /S^ (^a' + ah (n^r^ + (^^ ))^) + 6^ 
+2aP (a + lb (r^^ + (^r^^)) \ph) (27) 



Now, because \eh) is defined as having no component 
perpendicular to the n-fold axis, we can neglect all 
components of Rn"^ along that axis and replace each 
(^i'"^^ + (^i'^^y) term with 2cos (^) /(^^). Then, 
by applying both (22) and (23), we are left with 



Finally, for a plane wave we note that 
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IV. EXAMPLES OF LIGHT SCATTERING 



and the expression then becomes 

{e\p) = |EO|'(e,K) 



(28) 



The remaining two inner products in (27), 
{p\ (q^^^^^) \p) and {p\p), can be thought of as analo- 
gous cases of each other where /^^^^ and vice 
versa. So we wih begin by expressing {p\ (a^^^^) \p) 

in terms of {ph\ and \ph) = {M^^^^y^eoa^^^^ \eh) using 
(25) 



(p|(a(^^))-^b) = (p^|(a(6^))-^[a2+2a/36+/32a2+/3^fe^ 
+ (a/36+/32a6)(^(f^)n^6iv) 

+ (^L'^^n(f^O^)](A^^'''^)"'eoa(^^)|e^). 

Every term wih commute with (M^^^^) ^ (a 

result of (17)) and thereafter the permutation matrices 
can be absorbed by \eh) 

(p|(a(^^))-^b) = (p.|(c.(^^))-^(M(^^))-^[c.^ 

+2a/3(a+i6(^(f^) + (^(f^))^))]eoa(^^)|e^). 

Here the terms in the square parentheses are the same as 
in (27), and we have already shown that they reduce to 

l^p. As such {p\ (q^^^^^) \p) and, analogously, {p\p), 
can be expressed as 



{P\ 



{p\p) = mT{Ph\Ph)- 



(29) 



Substituting (28) and (29) into our original expressions 
for the extinction, absorption and scattering cross sec- 
tions (27) we get 



k 



lm{{eh\ph)} 



(30) 



{Ph\Ph) +Im 



\Ph 



Noticeably the cross sections are independent of the coef- 
ficients a and (3, which defined the generality of the inci- 
dent field in (24). As such, we can conclude that the ex- 
tinction, absorption and scattering cross sections will all 
be independent of the polarization in a linearly-polarized 
plane wave that is travelling parallel to the n-fold axis of 
a dipole system (for n > 3). 



To demonstrate the validity of our approach we em- 
ployed a number of methods to study the light scattering 
by some oligomer structures with n-fold symmetries [see 
Figs. 2-4]. First, we applied analytical Generalized Mie 
theory [38, 39] for spherical particle clusters to calcu- 
late the exact total extinction, scattering and absorption 
cross-sections. Second, we used CST Microwave Studio 
to calculate the near-field profiles of corresponding struc- 
tures. And, finally, we employed the dipole approxima- 
tion and dyadic Green's function method to analyze op- 
tically induced electric and magnetic dipoles distribution 
of individual nanoparticles [37]. We should mention that 
all three approaches prove the polarization independence 
of optical properties, in accordance with our theoretical 
predictions. Moreover, these methods provide us with 
flexibility in analyzing the structures while maintaining 
a full understanding of the system. We used Palik's book 
for permittivity data of various materials used in this pa- 
per [40]. 

Figure 2 show the extinction, scattering and absorption 
cross-sections for a trimer consisting of three touching sil- 
icon nanoparticles of radii R = 75nm. Figure 3 shows the 
extinction, scattering and absorption cross-sections for a 
heptamer consisting of seven gold nanoparticles with 6- 
fold symmetry. All particles of these oligomer-like struc- 
tures are in the same transverse plane, so the excitation 
field is identical for all particles. We can also lower the 
symmetry of a structure by shifting some particles along 
the propagation axis. For example, in Fig. 4 we present a 
structure with 3-fold symmetry which was derived from 
a gold heptamer structure shown in Fig. 3. To construct 
the new structure we, first, shifted two equilateral trimers 
of the outer ring in opposite directions from the central 
particle evenly spaced along the propagation axis, and, 
then, twisted with respect to each other. The final struc- 
ture appears to be chiral with 3-fold symmetry. Thus, 
according to our theoretical prediction we expect that 
it should be polarization invariant, which is numerically 
confirmed. 

All structures exhibit polarization-independent optical 
properties for any incident polarization angle, which does 
not necessarily coincide with the rotational symmetry of 
the structures. These figures show that the total absorp- 
tion is polarization-independent even though the near- 
field distribution varies with the incident polarization. It 
allows us to conclude, that, although, all structure ex- 
hibit some degree of geometrical anisotropy, their optical 
response is isotropic one. And the only requirement that 
we imply is that the structure supports n-fold symmetry 
with n > 3. 

Based on the coupled dipole approximation 
method [41] the optical response of structures with 
arbitrary geometries and complex refractive index can 
be approximated by an ensemble of discrete dipoles. 
Thus, our results can be easily generalized to any 
structure with n-fold symmetry. Figure 5 shows the 
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results of direct numerical simulations of a continuous 
structure with 3-fold symmetry for simplicity modeled 
as p{0) = R[l + cos(1.56>)] with R = 200 nm (on the 
transverse plane) and h = 100 nm (along longitu- 
dinal direction) thick made of gold. It still exhibits 
polarization-independent optical response, in full agree- 
ment with our approach above. This proves that our 
statement is quite universal and can be applied to any 
system. 



all rotational symmetry of a structure without any con- 
ditions placed on other ingredients of a given system. 
Specific numerical examples have been demonstrated 
to follow these derivations. Such robust polarization- 
independent features are expected to play an important 
role in various applications including nanoantennas, sens- 
ing, imaging, solar cells and many other applications in 
biological and medical research. 



CONCLUSIONS 
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FIG. 1. Schematic view of the problem under consideration - 
linear polarised plane wave scattering by a system with rota- 
tional symmetry. 
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FIG. 2. Extinction, scattering and absorption cross-sections 
for two orthogonal polarizations for trimer structure made of 
touching silicon nanopartices of radii R = 75nm calculated 
by using (a) discrete dipole approximation (DDA) method 
and (b) by direct numerical simulations with CST Microwave 
Studio. Bottom panels show magnetic field distribution at 
the Fano resonance [42], A = 612nm, for two polarizations 
together with induced magnetic dipole moments. 
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FIG. 3. Extinction, scattering and absorption cross-sections 
for two orthogonal polarizations for heptamer structure made 
of gold nanopartices of central Rc = 65nm and outer Ro = 
75nm radii separated by d = lOnm calculated by using (a) 
discrete dipole approximation (DDA) method and (b) by di- 
rect numerical simulations with CST Microwave Studio.. Bot- 
tom panels show electric field distribution at the Fano reso- 
nance [18], A = 600nm, for two polarizations together with 
induced electric dipole moments. 
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FIG. 4. Extinction, scattering and absorption cross-sections 
for a structure with 3-fold symmetry consisting of two equi- 
lateral trimers (marked by different colors) that are separated 
by a central particle and evenly spaced along the propagation 
axis. The trimers are also twisted with respect to each other 
in order to make the structure chiral. 
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FIG. 5. Scattering and absorption cross-sections for continu- 
ous structure for two orthogonal polarizations with n = 3-fold 
symmetry. 



